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Electron transport through multi-terminal rectangular arrays of quantum rings is studied in the 
presence of Rashba-type spin-orbit interaction (SOI) and of a perpendicular magnetic field. Using 
the analytic expressions for the transmission and reflection coefficients for single rings we obtain the 
conductance through such arrays as a function of the SOI strength, the magnetic flux, and of the 
wave vector k of the incident electron. Due to destructive or constructive spin interferences caused 
by the SOI, the array can be totally opaque for certain ranges of k, while there are parameter values 
where it is completely transparent. Spin resolved transmission probabilities show nontrivial spin 
transformations at the outputs of the arrays. When point-like random scattering centers are placed 
between the rings, the Aharonov-Bohm peaks split, and an oscillatory behavior of the conductance 
emerges as a function of the SOI strength. 
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I. INTRODUCTION 

Magnetoconductance oscillations of quantum rings 
made of semiconducting materials^ exhibiting Rashba- 
type spin-orbit interaction^ 3 -^ (SOI) have been in- 
tensely studied in the past few years. These effects 
are manifestations of flux- and spin-dependent quan- 
tum interference phenomena. In view of the possible 
spintronic applications and the conceptual importance 
of these interference effects in multiply-connected do- 
mains, closed single quantum rings (without attached 
leads) j5i§jL2i as well as two- or three-terminal ones were 
mvestigateM^^^^^^^^^ 2 ^^ extensively. 
Additionally, the conductance properties of a linear chain 
of rings have also been determined. 23 

In this paper we present a method that enables one to 
calculate the conductance and the spin transport prop- 
erties of two-dimensional rectangular arrays consisting of 
quantum rings with Rashba-type SOI 2 ^ and with a per- 
pendicular magnetic field. Such arrays, fabricated from 
e.g. an InAlAs/InGaAs based 2DEG^ have been stud- 
ied in a recent experiment^ and in a subsequent theoret- 
ical work>21 to demonstrate the time-reversal Aharonov- 
Casher effect^ Here we present a more general survey of 
the magnetoconductance properties of such devices, in- 
cluding the perturbative treatment of the magnetic field 
which still allows us to analytically solve the scattering 
problem in case of two-, three- and four-terminal rings, 
which are then used as building blocks of larger arrays. 
Our method is based on analytic results, and can be 
used for an arbitrary configuration. For the sake of defi- 
niteness, we consider 3x3, 4x4 and 5x5 rectangular ar- 
rays ; 26 ' 27 which are closed in the vertical, and open in the 
horizontal direction. Additionally, we study the magne- 
toconductance properties and spin resolved transmission 
probabilities of the same array geometry with only one 
input channel. We also investigate to what extent the 
conductance properties are modified by the presence of 



point-like random scattering centers between the rings. 
In our calculations we assume that the rings are narrow 
enough to be considered one-dimensional and the trans- 
port of the electrons through the arrays is ballistic. We 
determine the magnetoconductance in the framework of 
the Landauer-Biittiker formalism^ 

Rectangular array s 26 ' 27 -depending on the number of 
input leads- consist of two-, three- and four-terminal 
rings (see Fig. [I}, where the two- and three-terminal 
ones are situated on the boundary of the arrays as shown 
in Fig. [2] with or without the input leads displayed by 
dashed lines. The transmission and reflection properties 
of two- and three-terminal rings have been determined in 
previous worksi"^^ 3 ^^^^ but the effect of 
the magnetic field on the spin degree of freedom has not 
been taken into account for an arbitrary geometry. Ad- 
ditionally, the most general boundary condition that is 
required by this two-dimensional problem has not been 
investigated so far. Therefore in Sec. |TT] we first con- 
sider a perpendicular magnetic field as a weak perturba- 
tion, then, in order to account for all possible reflections 
and transmissions when building up the array from sin- 
gle rings, we generalize our previous results to the case 
when electrons can enter/exit on any of the terminals of 
a three-terminal ring (results for two- and four-terminal 
rings are presented in the Appendix). Next, in Sec. IIII Al 
the individual rings are used as building blocks of the 
arrays by fitting the wave functions and their derivatives 
in the points where neighboring rings touch each other. 
Magnetoconductance properties are presented here as a 
function of the wave number k of the incoming electron, 
the magnetic flux and the SOI strength. Spin resolved 
transmission probabilities on the output side of the arrays 
are also derived. In Sec. IIII Bl we investigate the effect of 
random Dirac-delta scattering potentials in between the 
rings. 
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II. BUILDING BLOCKS OF 
TWO-DIMENSIONAL ARRAYS: SINGLE 
QUANTUM RINGS 

In this section we consider a single, narrow quantum 
ring^i of radius a located in the xy plane in the presence 
of Rashba SOI 24 and a perpendicular magnetic field B. 
If B is relatively weak, then the interaction between the 
electron spin and the field, i.e. the Zccman term can be 
treated as a perturbation and the relevant dimensionless 
Hamiltonian read s 11 ! 36 



H = 



.d_ 



<^so 
2fl ' 



^so 

An 2 



pi 



(1) 



where tp is the azimuthal angle of a point on the ring, <!> 
denotes the magnetic flux encircled by the ring, <& = /i/e 
is the unit flux, and cu^o = oe/Ka is the frequency asso- 
ciated with the spin-orbit interaction, hfl = h 2 /2m*a 2 
characterizes the kinetic energy with m* being the effec- 
tive mass of the electron, and the radial spin operator is 
given by a r = a x cos ip + a y sin (p. The perturbative term 
Hp is given byii 



of the incident electron, which is described as a plane 
wave. Within this approximation, the eigenspinors are 
not perturbed and their direction is still specified by the 
angle 9 given by Eq. ([3]). The energy eigenvalues includ- 
ing the first-order corrections are given by 

-1 1 



w 



Imposing the condition of energy conservation k 2 a 2 
E^\k) determines the possible values of k: 



= (-1 



where fi.j ~ 1, 2 and 



q 2 + (-!)' 



Q w 



(4) 



with q = y/(cj so /2n) 2 + E/hfl, and E = h 2 k 2 /2m* de- 
noting the energy of the incoming electron. The corre- 
sponding four eigenspinors read 



e- lv / 2 cos{9/2)\ 
e^/ 2 sin(6</2) ) 



(5) 



H ? = IT**' 

where u>l = g*eB /Am with g* and m being the effective 
gyromagnetic ratio and the free electron mass, respec- 
tively. 

The energy eigenvalues of the unperturbed Hamilto- 
nian are 

E (%) = (« - ^H^f (« - »+\, (A* = M) , 

(2a) 

and the corresponding eigenvectors in the | | z ), | [ z ) 
eigenbasis of a z read 

i \ - fe~' Ltp/2 u { -^\ 

^ W) = (% /2 J M) J , ( 2b ) 

where w (1) = -u (2) =cos(<9/2); = u (1 ) = sin(0/2) and 

n 



tan(0/2) = (l-w), 



(3) 



with w = \J\+ ojgo/Q 2 . 

The matrix elements of H p in the basis of these eigen- 
statcs arc obtained as 



(m) 



vM = (-if +i ^ cose = ^- 



(1) 



In the first-order approximation one neglects the off- 
diagonal elements; this is reasonable if they are small, 
i.e., if lu-l/SI <C k 2 a 2 , where k denotes the wave number 



(2)/ (2) \ 



e -i<p/2 s i n (gi/2)\ 
-e^/ 2 cos(6>/2)y 
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The wave functions belonging to the same energy in the 
different sections of the ring are linear combinations of 
these eigenspinors. 

The building blocks of the rectangular arrays we in- 
vestigate are two-, three- and four-terminal quantum 
rings (see Fig. [IJ, where, in general, the boundary con- 
ditions allow both incoming and outgoing spinor valued 
wave functions at each terminal: 'J,; = fie tkXi + rie~ %kXi 
(i = I, II, III, IV), where Xi denotes the local coordi- 
nate in terminal i. Note that the amplitudes fi,rj,fjj, 

'(//) T N 



etc. refer to two-component spinors, e.g. fi = 



For the sake of definitencss, we focus on a general three- 
terminal ring, shown in Fig. 1 (b) The scattering problem 



in the case of a ring with four terminals (Fig. 1 (c) ) can 
also be solved analytically, as presented in the Appendix, 
where we also give the results for a general two-terminal 
ring (Fig. 1(a)). The outgoing spinors (ri, i = I, II, III) 
are connected to the incoming ones (fi) by 2 x 2 matri- 
ces, which can be determined by requiring the continuity 
of the wave functions, and vanishing net spin current 
densities (Griffith conditions) 11 ! 13 ' 32 ! 37 at the junctions. 
For the same boundary conditions as in Ref. [35|, i.e. for 
fn-i fin = in Fig. 1(b) the reflection matrix which con- 



nects 77 to the incoming spinor // is given by 



TT 

R 



n ~ 

r>fl _ 
nfl _ 



cos 2 (6>/2) + sm 2 (d/2) - 1, 
- g^)sm{0/2)cos{0/2), 

r>fl 

■"TJ.' 

Q W sin 2 (6»/2) + cos 2 ((9/2) - 1, 



(7) 
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FIG. 1: The notations used for the spinor part of the wave 



funct ions in the case of two- (a) three- (b) and four-terminal 
rings (c) 



where 



qM = 8ka/y^ {-i(^) 2 sin(2 g ^7r) 



-kaqM sin(9^7i) sm(q^ (2tt-ti)) 
+ sm(q^j 2 ) wa.(qM (2ir- l2 )) 
+ik 2 a 2 sin(q M 71 ) sin(g (Al) (72 - 71 )) 
xsin(gM(27r-7 2 ))}, 



and 
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-12ikc(q M ) 2 sm(2q M TT) + Ak 2 a 2 q (tl) cos(2q itl \) 
-2k 2 a 2 q W [cos(2g (Al) (77-72+71)) -cos^g^V) 

+cos(2g (Al) (7r-7 2 )) + cos(2g (Al) (7r-7i)) 

sin(2g (M) (tt - 72 +71 ) ) - sm(2q ( + l) tt) 

+sin(2g^)(7r-7i))-sin(2^)(7r-7 2 ))l , 



with <fi = Q/&o- The matrices describing the connection 
between the outgoing spinors rjj , Tin and the input // - 



the so called transmission matrices - are given by 



(^n ) || 
(^n0|| 
(-^n J )j.| 



e -W3 T (D 
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•^W^+r^sm 2 ^)), 
(r^-r^) sin(0/2)cos(#/2), 
3 <7n/a ( r W_ T W) sm (0/2) cos(0/2), (8) 
(rWsin 2 (0/2) + r( 2 )cos 2 (0/2)), 



where n = 1,2, indicating the two possible output chan- 
nels and 



= gW^ p i7i/2((-l)" +1 W +2^) 



{-fcasin( g ( ^(72 - 7 i))sin( 9 ( ' j) (2^ - 72 )) 

e -^((-l)" +1 -+20) sin(g (M) 7l) 

(^)(27r-7i))]}, 



— sin 
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{fca e - M ( ( - 1)M+1 ^)sin( g ^ 7l ) sin(qW( 72 - 7 i)) 
5in(^)(27r- 72 ))]}. 



- sin 



Note that the boundary conditions applied to obtain the 
R? 1 and T^ 1 matrices above are similar to that of Ref . l35l . 
However the magnetic field induced shift of the spin Zee- 
man levels leads to a doubling of the parameters accord- 
ing to Eq. ([4]). This modifies significantly the physical 
transport properties of the device. 

Let us point out that having obtained the matrix ele- 
ments above is enough to handle the problem with both 
incoming and outgoing waves on all terminals of the ring 
as shown in Fig. 1(b) Namely, we can consider the three 
inputs fi (i = I, II, III) separately and determine the 
corresponding reflection and transmission matrices. The 
outputs in the superposed problem will consist of contri- 
butions from all inputs: the reflected part of the spinor 
which enters on the same lead, and the transmitted parts 
of the other two inputs into the respective lead: 



ri = R fl fi + Ti n f n + T^'f m , 
rin = T^fx + T^fn + Rtefin. 



(9) 



Considering /// (fin) as the only input, the reflection 
and transmission matrices are the same as those for the 
input //, except for the appropriate changes of the an- 
gles, since in the reference frame of fn (fin), the angles 
of the output leads are measured from the lead through 
which fjj (fni) enters the ring. In order to get the con- 
tributions to the output spinors for the input /// (fin) m 
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the reference frame of //, the matrices need to be rotated 



(see Fig. 1(b) ) by the angle of 71(72) 



U' 1 

72-71 71 

72^27T-7i 



M f '" = U J2 M\^_ J2 U- 1 



(10) 



(11) 



72^27T-72+7l 



where M = R,T 1: T 2 and 









e 2 



n = 1,2. 



The above approach is also valid in the case of the two- 
and four-terminal rings. Using the reflection and trans- 
mission matrices as presented in the Appendix, the more 
general problem of having both incoming and outgoing 
waves on all terminals can easily be treated. All possible 
reflections and transmissions can thus be taken into ac- 
count when forming two-dimensional arrays of such rings. 



III. RECTANGULAR ARRAYS OF QUANTUM 
RINGS 



A. Magnetoconductance properties 

Based on the analytic results presented in the previ- 
ous section and in the Appendix we may build NxM 
two-dimensional rectangular arrays of quantum rings, 
where both perpendicular electric and magnetic fields are 
present, so that the former one can be used to change the 
strength of the SOL— Here we focus on of 3x3, 4x4, and 
5x5 arrays and assume that neighboring rings touch each 
other. In addition, we limit ourselves to arrays that are 
closed in the vertical, and open in the horizontal direc- 
tion, as shown in Fig. [5] Two types of such arrays will 
be investigated: i) the electron can enter/exit the array 
through any of the rings in the horizontal direction, ii) 
the electron can enter the array through one ring only 
(no leads are attached to the other rings on the entrance 
side) , but can exit through any of the rings on the oppo- 
site side (Fig. [2] without the dashed curves). In both cases 
the conductance is derived from the linear set of equa- 
tions resulting from the fit of the wave functions 'J,- 
(i = I, II, III, IV and k,l = 1..N, where N is the num- 
ber of rings along one direction in the array) and their 

derivatives d 00 in the points, where the rings touch 



each other, i.e. for example: 
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Here we used the notations of Fig. O (Note that the 
negative sign in Eq. (|12p is a consequence of the opposite 



direction of the local coordinates in the leads III of ring 
{11} and I of ring {12}) Eqs. JI2J) lead to 
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in 
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= ~fi 
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from which follows that 

f (n) 
J in 

r (n) 
'ill 



.(12) 



fi 



(12) 



i.e., the spinor entering (exiting) ring {11} on terminal 
III is equal to the spinor exiting (entering) ring {12} on 
terminal /. The spinors r\jP and r^ 12 -* can be given with 
the help of the reflection and transmission matrices of a 
three-terminal ring according to Eqs. 

For a small number of rings the resulting set of equa- 
tions can be solved analytically, however already for an 
array of 3x3 rings shown in Fig. [21 it consists of 60 equa- 
tions, which is preferably solved by numerical means, al- 
though analytic solutions exist in principle. (For larger 
arrays the number of equations scales practically with 
the number of rings). After having determined the out- 

JNN) 
11 ' 



put spinor valued wave functions , r , 
where N is the number of rings in the horizontal direc- 
tion, the Landauer-Biittiker— formula 



where 
G T = 

G.i = 



G — Gj + 



(r<r>) T 
(r<r>)i 



is used to calculate the conductance of the arrays, aver- 
aged over the two a z eigenspinor inputs. We note that 
our method of using single rings as building blocks can 
easily be used to determine the conductance of arrays of 
arbitrary - not necessarily rectangular - configuration as 
well. 

Figure [3] shows a contour plot of the conductance (in 
e 2 /h units) of rectangular arrays of 3x3, 4x4 and 5x5 
quantum rings, for zero magnetic flux as a function of the 
SOI strength wso/^ a- n d ka. The values of ka are var- 
ied around kpa — 20.4, corresponding to a Fermi energy 
11.13 meV in case of an effective mass m* = 0.023m of 
InAs and rings of radius a = 0.25 fiia. In two-dimensional 
electron systems within an InAs quantum well, the value 
of a can be varied^ up to 40 peVm. The different ar- 
rays show similar behavior for larger values of the SOI 
strength: there are slightly downwards bending stripes 
(initially around even values of ka) where the devices are 
completely opaque for the electrons, and also conducting 
regions which are initially around odd values of ka and 
have complex internal structure. Comparing our results 
to the case of a single ring with diametrically coupled 
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FIG. 2: The geometry of the device in the simplest case of 
a 3x3 array with three or one (without leads displayed with 
dashed lines) input terminals. The notations can easily be 
generalized to larger arrays. 



leadsii, it can be seen that the overall periodicity as a 
function of ka is determined by single-ring interferences. 
The increasing number of the rings causes modulations 
superimposing on the single-ring behavior. This point is 
probably the most apparent if we recall^ that zero con- 
ductance areas are simply lines on the ka — u^o/Vl plane 
for a single two-terminal ring, while in our case there are 
stripes, the width of which is slightly increasing with the 
size of the array. This effect is related to the increasing 
number of consecutive partially destructive interferences 
that finally lead to essentially zero currents at the out- 
puts. Additionally, if we considered an infinite network, 
the periodic boundary conditions would allow only dis- 
crete values of ka for a given SOI strength with nonzero 
conductance. Thus the results presented in Fig. [^demon- 
strate a transition between the conductance properties of 
a single ring and that of an infinite network. 

Focusing on small values of ui^o/^l, Fig.[3]shows a nar- 
rowing of the non-conducting regions until they eventu- 
ally disappear when no SOI is present. Here the conduc- 
tance still depends on ka, but its minimal values are not 
zeros and a periodic behavior can be seen: for a network 
of N x N rings, there are N minima as the value of ka is 
increased by 1. This size-dependent modulation is related 
to the horizontal extent of the device: If we compare the 
conductance of the networks to that of rings of the same 
size and number without vertical connections, the same 
periodic behavior can be seen around zero SOI. 

Figure U shows the normalized magnetoconductance 
of networks of 3x3, 4x4 and 5x5 quantum rings for 
ka = 19.6 as a function of the SOI strength and the mag- 
netic flux <E> (measured in units of <!>o). When wso/^ 
is zero, Aharonov-Bohm (AB) oscillations appear. For 
larger values of u>so/^ both AB and Aharonov-Casher— 
oscillations can be seen in the magnetoconductance. As 
Fig. 0] was plotted for a certain value of ka, the effect 
of the bending non-conducting stripes seen in Fig. [3] can 
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FIG. 3 : (c olor o nlin e) The co ndu ctance G/Go (Go = e 2 /h) of 
a 3x3 (a) 4x4 (b) and 5x5 (c) rectangular array with 3, 4, 
and 5 input terminals, respectively, for zero magnetic flux as 
a function of the SOI strength and ka. 



also be seen as the decrease of the conductance when 
such a stripe is reached due to the change of the SOI 
strength, and its increase again, when the stripe is left. 
We note that for larger values of ka this bending effect 
is less pronounced. 

Figures [5] and [6] show the conductance of a 5x5 net- 
work with a single input lead in the middle (i.e. attached 
to ring {31}, using the notations of Fig. [5]) as a function 
of ka, and tuso/Q (Fig- E]), an d the magnetic field and 
^so / fi (Fig. [6]) . The overall structure of these plots re- 
mains the same as in the case when the current can enter 
through all the rings on the left hand side, but the dif- 
ferent boundary conditions modify the fine structure of 
the plots. 

Our method allows the calculation of the spin direc- 
tions for the different output terminals, and we found 
that spin-dependent interference in the array results in 
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FIG. 4: ( colo r onli ne) T he cond ucta nce G/Go (Go = e 2 /h) 



of a 3x3 [(a)] 4x4 
3, 4, and 5 input terminals 



(b) and 5x5 (c) rectangular array with 
respectively, for ka 



19.6 as a 

function of the SOI strength and the magnetic flux $ (in units 
of 4>o = h/e). 



nontrivial spin transformations. Fig. [7] shows the spin 
resolved transmission probabilities for a 5x5 ring array 
with a single input lead. The incoming spin state is cho- 
sen to be | 1 Z ), i.e., the spin-up eigenstate of a z , and the 
contour plots show the probabilities of the | 1 X ) 7 | f y ) 
and | tz) outputs at ring {55} on the right hand side. 
The fact that the | | z ) input spinor changes its direction 
(as it is seen in Fig. it can be transformed into | tz) or 
1 y )) is due to the SOI induced spin rotations. The ac- 
tual values of the spin resolved transmission probabilities 
are determined by the spin dependent interference phe- 
nomena. Fig.[5]shows the z component of the normalized 
output spinors and visualizes that spin resolved results 
depend on the input side geometry as well. As we can 
see, the spin components change in the whole available 
range between -1 and 1, and their behavior is rather dif- 
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FIG. 5: (color online) The conductance G/Go (Go = e 2 /h) 
of a 5x5 rectangular array with a single input lead attached 
to ring {31}, for zero magnetic flux as a function of the SOI 
strength and ka. 
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FIG. 6: (color online) The conductance G/Go (Go = e 2 /h) 
of a 5x5 rectangular array with a single input lead attached 
to ring {31}, for ka = 19.57 as a function of the SOI strength 
and the magnetic flux $ (in units of $o = h/e). 



ferent for the cases when the electron can enter the array 
through any of the five terminals, or only through the one 
attached to ring {31}. This phenomenon together with 
other spin dependent interference effect o 38 ' 39 i 4[l ' 41 i 42 ' 43 i 44 
can lead to spin sensitive quantum networks. 



B. Effect of point-like scatterers 

Now we will investigate to what extent the conduc- 
tance properties are modified by the presence of ran- 
dom scatterers. Although high mobility samples have 
already become available (such that at cryogenic tem- 
peratures transport is found to be ballistic over tens of 
microns), considering also the effects caused by scatter- 
ing events provides a more realistic description for most 
cases. To this end we introduce point-like scattering cen- 
ters between the rings. In particular, at each point j 
where two rings touch each other, we consider an ad- 
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FIG. 7: (color online) The probabilities of the | 1 X ) (a) | j^) 



(b) and | | z ) (c) outputs at ring {55} of a 5x5 rectangular 
array with one input lead (attached to ring {31}), for ka = 
19.6 as a function of the SOI strength and the magnetic flux 
$ (in units of $0 = h/e). The incoming spin state is chosen 
to be I Ta). 



ditional Dirac delta potential of the form rjj5(j). Here 
rjj represent independent normally distributed random 
variables with zero mean, and root-mean-square devia- 
tion D. By tuning D we can model weak disturbances 
(small D) as well as the case when frequent scattering 
events completely change the character of the transport 
process (corresponding to large values of D). 

As shown in Fig. [9j the most general consequence of 
these random scattering events is the overall decrease 
of the conductance. However, for strong enough dis- 
turbance, more interesting effects can be seen, namely 
the splitting of the AB peaks which is more pronounced 
with increasing D. Note that the scattering has the most 
dramatic effect for the AB resonances, i.e. $ = n$o 
and the least for the anti-resonance condition, i.e. $ = 
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FIG. 8: The spin transformation properties of a 5x5 array 
with input leads attached to all rings and only to ring {31} 
(black and grey curves, respectively). The z component of the 
normalized spin states transmitted via the output terminals 
attached to ring {25} (solid line) and ring {45} (dashed line). 
The incoming spin state is chosen to be j | z ). 
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FIG. 9: (color online) The conductance G (in units of Go = 
e 2 IK) of a 5x5 rectangular array with point-like random scat- 
terers between the rings for different root-mean-square devi- 
ations D as a function of the magnetic flux $ (in units of 
$0 = h/e) for ka = 20.2 and uj S o/M = 13.0. 



(n + 1/2) <f> - Fig. HOI shows that the conductance has an 
oscillatory dependence on the strength of the SOI. No- 
tice that the introduction of scattering reduces the pe- 
riodicity with a factor of two. We want to stress that 
the model we considered (random elastic scattering pro- 
cesses in single-electron approximation) is similar to the 
case when the Altshuler-Aronov-Spivak (AAS) effect^ 
is expected to survive in a single ring. Our results for 
a more complex geometry indicate similar physical con- 
sequences of the scattering events: introduction of new 
peaks in the AB oscillations and appearance of u>so/^ 
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FIG. 10: (color online) The conductance G (in units of 
Go = e 2 /^) of a 5x5 rectangular array with point- like random 
scatterers between the rings for different root-mean-square de- 
viations D as a function of the SOI strength for ka = 19.6 
and $ = 0.3$ . 



Appendix 

Here we present the detailed analytic expressions of the 
scattering problem for a general two- and four-terminal 
ring, in which SOI and a perpendicular magnetic field 
is present, the latter of which is considered as a pertur- 
bation. As we have shown in Sec. [TTJ it is sufficient to 
consider only one input terminal and determine the con- 
nection between the input and output states, i.e. the 
reflection and transmission matrices, since the more gen- 
eral boundary condition of having inputs on all terminals 
is just a superposition of such cases with an appropriate 
rotation of the matrices (see Eqs. (fT0^ - (lTTj) ). Consider- 
ing fx as the only input (i.e. f^x — 0, in Figs. 1(a) and 
l(c)| , requiring the continuity of the wave functions and 
applying Griffith boundary condition a 32 ! 37 at the junc- 
tions in both cases, we can obtain the reflection matrices 
W* and W 1 of the two-terminal ring, and of the four- 
terminal ring, respectively. Both can be written in a form 
analogous to that of Rf* of the three terminal case given 
by Eq. (0) with 



conductance oscillations for a given value of 4>. This lat- 
ter effect has recently been predicted for a single ring^ 
and was detected in the case of ring arrays^ 



4k 2 a 2 
+iq M sin 



{sin(g ( ^ 7l ) sin( ? W(27T-7i)) 
in(2 g ( ^7r)| , 



IV. SUMMARY 



and 



In this paper we calculated the spin dependent trans- 
port properties of two-dimensional ring arrays. We ap- 
plied general boundary conditions for the case of single 
quantum rings, which allowed the construction of arrays 
of such rings as building blocks. The magnetoconduc- 
tance of two-dimensional arrays of 3x3, 4x4, and 5x5 
quantum rings exhibited Aharonov-Bohm and Aharonov- 
Casher oscillations. We also determined the spin resolved 
transmission probabilities of the arrays and found signif- 
icant spin rotations depending on the SOI strength. We 
introduced point-like random scattering centers between 
the rings, which, for strong enough disturbance resulted 
in the splitting of the AB peaks and the emergence of an 
oscillatory behavior of the conductance, characteristic to 
the SOI strength. 
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2ka r 
yM X 



fc 3 a 3 



cos(2gHr) + cos(2q (Al) (7r-7 3 + 72 - 7i)) 

-cos(2g^ ) (7r- 73 + 72 )) + cos(2 (? ( ' j) (7r-73+7i)) 

-cos(2 ? W(7r-7 2 +7i))-cos(2gW(7r-73)) 

+cos(2(7 (Al) (7r- 72 ))-cos(2q ( ' i) (7r- 7 i)) 

+2ik 2 a 2 q M [sin(2q (Al) (7r-73+7 2 ))-3sin(2<7^V) 

+sin(2g ( ' l) (7r-73+ 7 i))+sin(2g^ ) (7r-72+7i)) 
+Aik 2 a 2 q (ti) [sm(2<? ( ^7r-7i))-sm(2g ( ^7r-7 3 )) 
-4fca(q ( ' l) ) 2 [cos(2q^(7r-7 3 ))+cos(2 g ^ ) (7r-7 2 )) 
+cos(2g ( ^(7r-7i))-3cos(2g^V) 



-8*(g 



sin(2 9 (Al) 7r)} 



respectively. Here 

t/M = fcV [cos(2g^)(7r-7 1 ))-cos(2g^) 7 r) 
+4:ikaq {tl) sin(2g W 7r) 

-4(q^) 2 [cos[((-lf +1 w + 2(t>)n\ +cos(2gHr) 



3 (2gHr) 



i ((-1)" +1 u;+20)7r - 
~32iA:a((j (Al) ) 3 sin(2g ( ^7r) + 24fc 2 a 2 (g^ ) ) 2 cos(2q ( ^7r) 
^tfc 2 a 2 (g^ ) ) 2 [cos(2q ( ^(7r-73))+cos(2q ( ^(7r-7 2 )) 
+cos(2q ( ^(7r-7i))+cos(2g ( ^(7r-7 3 +7i)) 
+cos(2q ( ' l) (7r-73+7 2 ))+cos(2( ? (Al) (7r-72+7i)) 
-8ifcVg (/l) Sm(2gHr) 

+4i/c 3 a 3 q ( ^ ) [sin(2( ? ^ ) (7r-73 + 72))-sin(2g^ ) (7r-73)) 
+sin(2g ( ^(7r-72+7i)) + sin(2 9 ( ' l) ( 7 r-7 1 )) 



cos(2g ( ' t) (7r-73 + 72 - 7i))+cos(2q (Al) 7r) 

+cos(2< ? ( ' i) (tt -73+71)) -cos(2g (Al) (7T-73+7 2 )) 
-cos(2q ( ' l) (7r-72+7i))-cos(2g ( ^(7r-7 3 )) 
+cos(2q ( ' l) (7r-7 2 ))-cos(2g (Al) (7r-7i)) , 



where the angles 7^ are defined in Fig. 1(a) and (c) The 



transmission matrices T^' of the two-terminal ring and 
(n = 1,2,3) of the four-terminal ring, can be given 
in an analogous form to that of the transmission matrices 
Tl 1 of the three-terminal one given by Eq. © with 



}0) 



sin( g ( ' j) (27r-7 1 )) 



-,.((-1)^^+20) sin(g (M) 7l) 



and 



\ik 2 a 2 sin( g ( ^(27r-2 7 3 + 272-7i)) 
-sin( 9 ( ^(27r-7i))+sin( g ( ^(27r-27 2 -|-7i)) 
-sin(^(27r-2 7 3 + 7i)) 
-2kaq M [cos(<7^ ) (27r-272 + 7i)) 
-2cos((7 ( ^(27r-7i)) + cos(g (Al) (27r-273+7i)) 
+4i{q^f [ e -"((-D f ' +1 -+20) sin ( ? M 7l ) 

-sin( 9 ^(27r- 7 i))] } , 



» _ 



4fcagM ai72/2 ((_i)M+i w+2<A ) 



^((-1)^+^+20) cos((7 (M) 72 ) 

^-^((-DM+^+a*) cos(ffM(2 7l -72)) 
+cos(g ( ^(27r-273 + 72))-cos(g^(27r-7 2 )) 
+4z(^)) 2 [ e — ((-D" +1 -+20) sin((? („) 72 ) 

-3111(^(2^-72))] } , 



sin(q (/i) 73 ) 



+ sin( (Z ^ I \27 1 -73)) -sin(g^(272 -73)) 
+ sin(< ? ( ' l) (272-27i-7 3 )) 
-2kaq^e-™ ((-D^ 1 ^) [ 2 cos(?°S) 
- cos(g^)(2 7l -73)) -003(5(^(272 -73)) 
+4 1 (< Z ^)) 2 ((-D^ +1 -+20) sin ( e (M) 73 ) 

-sin(^)(27r-7 3 ))]}, 



respectively. 
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